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We propose a physically transparent analytic model of astrophysical S-factors as a function of a 
center-of-mass energy E of colliding nuclei (below and above the Coulomb barrier) for non-resonant 
fusion reactions. For any given reaction, the S(.E)-model contains four parameters [two of which 
approximate the barrier potential, U(r)]. They are easily interpolated along many reactions in- 
volving isotopes of the same elements; they give accurate practical expressions for S(E) with only 
several input parameters for many reactions. The model reproduces the suppression of S(E) at low 
energies (of astrophysical importance) due to the shape of the low-r wing of U (r). The model can be 
used to reconstruct U(r) from computed or measured S(E). For illustration, we parameterize our 
recent calculations of S(E) (using the Sao Paulo potential and the barrier penetration formalism) 
for 946 reactions involving stable and unstable isotopes of C, O, Ne, and Mg (with 9 parameters 
for all reactions involving many isotopes of the same elements, e.g., C+O). In addition, we analyze 
astrophysically important 12 C+ 12 C reaction, compare theoretical models with experimental data, 
and discuss the problem of interpolating reliably known S(E) values to low energies (E < 2 — 3 
MeV). 

PACS numbers: 25.70.Jj;26.50.+x;26.60.Gj,26.30.-k 



I. INTRODUCTION 

Nuclear reactions are very important [J-Q for the 
structure, evolution, nucleosynthesis and various obser- 
vational manifestations of main-sequence stars, giants 
and supergiants, presupernovae, white dwarfs and neu- 
tron stars. Depending on temperature, density and other 
parameters, stellar burning may involve many reactions 
of different nuclei, from light to heavy, and from stable 
to neutron- and proton-rich ones. Their rates can be cal- 
culated using the reaction cross sections a(E), or related 
astrophysical iS-factors defined as 

o(E) = E- 1 exp(-27rry) S{E). (1) 

Here, E is the center-of-mass energy of the rea ctants 
[(Ai.Zi) and {A 2 ,Z 2 )}, v = a /(hv) = y/E R /E is the 
Sommerfeld parameter, v = y / 2E/fi is the relative ve- 
locity of the reactants at large separations, a = Z\Z 2 e 1 , 
Er = a 2 /j>/(2K ) is similar to the Rydberg energy in 
atomic physics, and fi is the reduced mass. The fac- 
tor exp(— 2irr]) is proportional to the probability of pen- 
etration through the Coulomb barrier U(r) = a/r with 
zero angular orbital momentum, assuming that this pure 
Coulomb barrier extends to r — > (for point-like nu- 
clei); E~ x factorizes out the well-known pre-exponential 



low-energy dependence of a(E). The advantage of this 
approach is that S(E) is a much more slowly varying 
function of E than exp(— 2ini) and o~(E). 

For astrophysical applications, one needs to know S(E) 
for many reactions at low energies, E < a few MeV. Ex- 
perimental measurements of cr(E) at such energies are 
mainly not available (because the Coulomb barrier expo- 
nentially suppresses low-energy cross sections). Theoret- 
ical calculations are model dependent, so that nuclear- 
physics uncertainties of calculated S(E) can be substan- 
tial. Theoretical calculations show that S(E) can vary 
by several orders of magnitude in the energy range of as- 
trophysical importance for a given reaction, and it can 
vary over many orders of magnitude for different reac- 
tions (e.g., [J, [j| and references therein). It is the aim 
of this paper to propose (Sec. UTJ) a physically transpar- 
ent analytic model of S(E) for non-resonant reactions 
between heavy nuclei in order to explain these features, 
simplify the use of available S(l?)-data, and clarify the 
problem of interpolating reliably known S(E) values to 
low energies of astrophysical importance. 

Astrophysical S'-factors have been parameterized by 
different analytic formulae (see 0, HH?]], f° r references). 
We think that our new model is more flexible. It al- 
lows one to approximate S(E) for many reactions with 
minimum number of fit parameters. For instance, our 
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recent approximation of S(E) for 946 reactions involv- 
ing different isotopes of C, O, Ne, and Mg with 8514 fit 
parameters is now replaced (Sec. IHIj) with the approxi- 
mation containing 90 parameters. Moreover, our model 
directly relates S(E) with the parameters of the effective 
potential U (r) of nucleus-nucleus interaction, and helps 
to reconstruct (constrain) U(r) from the S(E) data (com- 
puted or experimental ones; see Sees. IHII and IIVI) . 



II. ANALYTIC MODEL 
A. General approach 

Let us construct a simplified model of S(E) at suf- 
ficiently low energies E at which the main contribu- 
tion to the reaction cross section comes from the s-wave 
channel. For the reactions involving the nuclei with 
Zi,Z2 ~ 6 — 12, this is true at the energies of a few 
tens of MeV. According to the theory of inelastic scat- 
tering (e.g., Ref. 0), a reaction cross-section at such low 
energies has different energy dependence above and be- 
low the Coulomb barrier. Below the barrier (at E < Ec, 
Ec being the barrier height) it behaves as 



<t{E) = S E- 1 exp$(E), 

HE) = -IjJdr^^U-E), 



(2) 
(3) 



where <$>(E) is the semi-classical exponent argument in 
the expression for the barrier penetrability. We adopt the 
semi-classical approximation to calculate this penetrabil- 
ity; U = U(r) is the effective nucleus-nucleus potential 
(it is Coulombic at large separations r but is affected by 
nuclear forces at small r); r% and r-i are classical turn- 
ing points. In Eq. (|2|), So is a slowly varying function of 
E which we treat as a constant. Its order-of-magnitude 
estimate in terms of physical quantities can be deduced, 
for instance, from the consideration in Appendix C of 
Ref. 0]: 



So 



2irh 2 



Ec 
V ' 



(4) 



where Vq ~ 40 MeV. The pre-exponent factor E~ x in 
© can be written as (1/VE)(1/VE), where one factor 
1 / \f~E~ is a generic feature of low-energy reaction cross 
sections (neglecting barrier penetration); the extra fac- 
tor \j\f~E comes from three-dimensional penetrability 
through the Coulomb barrier Let us stress that So 
is not identical to the astrophysical 5(S)-factor. We will 
see that So, in contrast to S(E), weakly depends on spe- 
cific reaction. 

At E > Ec the Coulomb barrier is transparent (in the 
semi-classical approximation), $(E) = 0, and a(E) oc 
1/yfE. 

Combining the definition of S(E), given by Eq. (fl}, 
with Eq. p]). and using the above arguments, we present 
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FIG. 1: (Color online) Model effective potential U(r) for the 
46 Mg+ 46 Mg reaction discussed in Sec. IIIII The solid curve is 
the pure Coulomb segment; the dot-and-dashed curve is the 
inverted parabolic segment; the dotted line shows the pure 
Coulomb potential extended to r — > 0. The filled dot marks 
the barrier hight (U — Ec at r = Rc), and the open dot 
separates the Coulomb and parabolic segments (U = Eci at 
r = Rci). 



S(E) in the form 

S(E) = So exp*(S), 

if>(E) = 2tti] + $>(E) at E < E c , (5) 
S(E) = So exp (27T7?) y/E/Ec 

x[l + C(E-E c )/E] at E>E C . (6) 

The last equation is phenomenological and contains a 
constant parameter £. This equation extends Eq. ([5]) to 
the energies above the barrier. At energies E ^> Ec 
(but still low enough for the model to be valid) we have 
S(E) = S exp(27T7 ? )(^ + l)y/E/E c . Thus, $ determines 
the magnitude of the reaction cross section at E 3> Ec; 
it is also important for describing the S(E) behavior at 
E ~ E C - 

B. Model of barrier potential 

Let us adopt a simple model of U(r), 
a 



Mr) 



U(r) = E c 



at r > Rci, 



1-/3 



(r - Rc) 2 



R c 



at r < Rci- (7) 



It is a pure Coulomb potential at r > Rci and an in- 
verse parabolic potential at smaller r. The parabolic 
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segment truncates the effective interaction at small sep- 
arations; Ec = U(Rc) is the maximum of U(r) (the 
barrier height). This model is the most natural and sim- 
ple approximation of U (r) which allows us to analytically 
calculate the barrier penetrability. We require U(r) and 
its derivative be continuous at r = Rci, and introduce 
<5 = (Rci — Rc)/Rc that characterizes the width of the 
peak maximum of U(r). 

In this way, U (r) is determined by two parameters, say, 
Ec and S, with 



Q (2 + 3<5) 
° 2£ c (l + <5) 2 ' 
R cx = Rc{l + 5), 

Eci — U(Rci) = Ec 



1 



.5(2 + 3(5)' 
2(5 



2 + 3(5' 



(8) 



The potential U(r) passes through zero at r = Rco = 
Rc (1 — /3 -1 / 2 ). Its behavior at smaller r does not affect 
directly our results. Realistic models should correspond 
to /3 >■ 1 (the low-r slope of U (r) should be sharp; Rco 
should be positive) which translates into 8 <C | (because 
f) = 1 corresponds to 8 = |). 

For example, in Fig.[T]we plot a model potential U(r) 
for the 46 Mg+ 46 Mg reaction. It will be discussed in Sec. 
Hm In this case, £ c =16.27 MeV and (5=0.0332, so that 
Eq. © yields i? c =12.53 fm and i? C i=12.95 fm. We show 
the pure Coulomb and parabolic segments (the solid and 
dash-dot lines, respectively) separated by the open dot. 
The filled dot is the potential maximum. The dotted line 
is the pure Coulomb potential extended to r —> 0. The 
thin solid and dashed vertical lines position, respectively, 
the maximum and separation points. Although 8 is for- 
mally small, it produces a noticeable U(r) wing at low 
r. 

With the potential the integral ([3]) is taken analyt- 
ically. At E < Ec\ we have 



Vr(E) 



* r (£) + (9) 
4 \ ~r ^arcsin ^/av + yf x r (l — x r )^j , 




E R {Ec - E) 
Ec Ec 



I — + arcsinx; + Xi\J 1 — xf j . 



where 7 = (2 + 3<5) 3 / 2 V^/(l + 8) 2 ; $ r (E) and $>i{E) 
contain the contributions from the integration regions 
of Rci < r < T2 and r± < r < Rci, respec- 
tiv ely; x r = Rci/ r 2 = RciE/a = E/E C i, and xi = 
8y/ ' j3Ec /(Ec — E). The term 2irq exactly canceled the 
opposite term which appeared after the integration in 
$(E). At E C i < E < E c we have 



*(J5) = 27r77 + *;(£') 




Ec — E 
Ec 




Equations © and ([10) fully determine *(£') in (JSJ in 
an analytic form. Then Eqs. ([5]) and ([S]) give an analytic, 
easily computable model for S(E). It contains four pa- 
rameters, So, Ec, 8 and £; each parameter has simple 
physical meaning. 

For astrophysical applications, one needs S(E) at sub- 
barrier energies. In this case, it is natural to present 
S(E) in the form 



S(E) = 5(0) eM9iE + g 2 E 2 + ...), 



(11) 



where g\, (/2,---are some expansion coefficients (deter- 
mined by low-r behavior of U{r); see below). The main 
energy dependence of S(E) at E < Ec is thought to be 
given by exp(<7i_E). To reduce the energy dependence one 
often introduces (e.g., @, 0|) the modified S'-factor 



S(E) = S(E)eM-9iE), 



(12) 



that is a much less variable function than S(E); the mod- 
ified 5-factor is usually treated as constant. Our model 
differs from the traditional approach: instead of S we 
prefer to introduce So- Both quantities are nearly con- 
stants (at E < Ec) but So changes within much narrower 

limits than S for different reactions (Sec. IIII[) . 

In our model, we can use Eq. ([5]) and expand the ex- 
ponent argument in powers of E. Keeping three lowest 
expansion terms, at E below Ec we come to Eq. (111[) 
with 



S(0) - S(0) = So exp( 5o ), 



(13) 



go being the zero-order expansion coefficient (followed 
by <7i, <72>- • • )■ The coefficients can be presented as <?.; = 
9ir + 9il (i = 0,1,2,.. .); gi T and gu collect, respectively, 
the contributions from the right (Coulombic, r > Rci) 
and left (parabolic, r < Rci) segments of the barrier 
potential. We obtain 




with xio = y/S/(2 + 3(5). We can write g = 
^Er/Ec MS), 9i = -y/ER/E^nid), 52 = 
y/ Er/E c (f2(8), where <fi(8) are functions of the only 
one argument 8. For instance, in the limit of 8 -C 1 
the first two functions are ^o(<5) = 8 — ir\/28 — 28 and 
ipi{8) = i-W2S-5. 
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The expansions like (fTTj) have been written long ago 
(e.g., Q and references therein) but only for the sharply 
truncated (rectangular - ret) Coulomb potential, 

U Ict {r) = - at r>R c , U Tct (r) = -V at r < R c . 
r 

This potential is obtained from our potential ([7]) in the 
limit of 8 — > 0, Eci — > Ec, and Rci — > Rc- In this case 
one usually used Eq. (jTTJ) neglecting c/2 and higher-order 
terms: 

S Ict (E) = 5 rct (0) exp(<? lrct £), 

fcct = *X = -^Vto ' (17) 

with S rct (0) = S rct = 5orct exp(g r0 ). Now £Jd = E c — 
a/ Rci; Rci should be treated as an effective radius of 
nucleus-nucleus interaction; the parameters gi r are not 
well defined because the radius Rc = -Rci of sharp bar- 
rier truncation is unphysical. One has <7i rc t < 0, that is 
Srct (E) always decreases with increasing E in the model 
(fT5| . We will show that it is more reasonable to use more 
advanced models of U(r) with the broadened peak and 
5i = 9\r + 9u- 



C. General properties of model S-factors 

The proposed model for S-factors is simple and can be 
analyzed in general form. 

a. Super-barrier energies. The model S- factor at 
E > Ec is given by Eq. ©. It is determined by three 
constants, Eq, So and £, being independent of 8. This 
S-factor rapidly decreases with the growth of E. 

Nevertheless, at E > Ec it is better to use the reaction 
cross-section instead of S(E). This cross section reads 
a(E) = S (E c E)- 1 ' 2 [1 + £(E - E c )/E]; it is a slowly 
varying function of energy. Let us stress again that our 
Eq. © cannot be extended to very large E, because the 
approximation of energy-independent So would become 
questionable. We expect that the term containing £ phe- 
nomenologically accounts for the contribution of higher 
partial waves t > at E > Ec- 

b. Subbarrier energies. In this case, S(E) is deter- 
mined by Sq, Ec and 8, being independent of £. Its main 
features are illustrated in Fig. [2j The left panel demon- 
strates the effective potential U(r). The vertical scale 
gives U(r)/Ec, which is independent of Ec- The hori- 
zontal scale shows t/Rq , where = a/ Ec for a given 
Eq- The solid line gives the sharply truncated Coulomb 
potential. The long-dashed, short-dashed, and dotted 
lines show U(r) for 8 =0.01, 0.1, and |, respectively. The 
higher-r wings of U(r) are pure Coulombic, U(r) = a/r, 
independent of 8. If Ec is fixed but 8 increases, the po- 
tential is broadened at small r and becomes, naturally, 
less transparent. The value 6 = ^ is critical, because 
nucleus-nucleus binding at small r disappears [we would 



have U(r) > at all r > 0]. This case is unrealistic for 
exothermic nuclear reactions (with positive Q-values) of 
our study. We show it as a limiting case. 



The middle panel of Fig. [2] presents 
^/Ec/Er \n[S(E)/S ] as a function of E/E c for 
subbarrier energies E < Ec [to visualize the shape of 
S(E)]. As seen from Eqs. (|9]) and (fTUj) . the presented 
function depends on the only one parameter, 8. The 
thick lines are calculated from Eqs. © and (Tit))) at the 
same values of 8=0, 0.1, 0.01, and ^ as in the left panel. 
Recall that for the sharply truncated Coulomb potential 
S(E) increases as E — > 0. For larger 8, the barrier is less 
transparent [which can lower S(E) to a great extent]. 
Moreover, the S(E) shape becomes different because the 
suppression of S(E) at lower E is naturally stronger. 
In the range of 0.05 < 8 < 0.15 the 5-factor turns 
into a slowly varying function of E, and at higher 8 it 
turns into a function which decreases noticeably with 
decreasing E, showing a pronounced hindrance of S(E) 
at low E. 



The thin lines in the middle panel of Fig. [2] present 
the same function, as the thick lines, but calculated us- 
ing the \ow-E asymptote (jlip . This asymptote is seen 
to be remarkable accurate for all values of 6 < | and 
for all energies E below Ec- Actually, the quadratic 
term <?2 E 2 in (fTTj) is small; it is often sufficient to use 
the reduced expression S(E) = So exp(<?o + g\E). The 
quadratic term modifies the exponent argument maxi- 
mum by 3% at E = E c and 5 = 0. Using Eqs. (QU and 
(JT5J) one can show that g\ changes sign at 8 = 0.07745 
(from minus at smaller 8 to plus at higher) . This explains 
a slow energy dependence of S at 0.05 < 8 < 0.15. At 
8 = \ the function In S(E) is almost linear in E, so that 
the thick and thin dotted lines nearly coincide. 



The right panel of Fig. [5] demonstrates log 10 [S , (0)/S , ] 
as a function of -JEr/Ec for the same four values of i5. 
For a fixed value of Er/Ec, the S(0)/So ratio depends 
only on 8, in our model. The displayed range of a/ Er/Ec 
corresponds to the reactions which we consider in Sec. 
IIIII The plot shows the values S(0)/Sq = exp(g ) which 
determine the zero-energy S-factor, S(0), important for 
astrophysical applications. For the pure Coulomb po- 
tential extended to r — > 0, we would have 5(0) = So- 
The potential cutoff at low r drastically enhances the 
barrier penetrability, and, hence, 5(0) [as well as, gen- 
erally, low-energy S(E)] over So, up to ~100 orders of 
magnitude! The highest enhancement is provided by the 
sharply truncated barrier (|16|) . If we fix Eq (and thus 
Er/Ec) an d increase 8, the enhancement can be strongly 
reduced. Some examples are given below. 
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FIG. 2: (Color online) ie/t: Normalized effective potential U(r)/Ec versus normalized radial coordinate. Mid: Normalized 
]n[S(E)/S ] versus E/E c [thick lines - full model, thin lines - Eq. |TT|]. Right: log 10 [S(0)/So] versus ^JE R /E C . Solid, 
long-dash, short-dash and dotted lines refer to 5—0, 0.01, 0.1, and 1/3, respectively. 



III. EXAMPLE: 946 REACTIONS INVOLVING 
C, O, NE AND MG ISOTOPES 

A. Calculations 

For example, let us consider a set of S-factors, which 
we have calculated recently [5| for fusion reactions in- 
volving various isotopes of C, O, Ne, and Mg. The cal- 
culations include stable, proton-rich, neutron-rich, and 
very neutron-rich isotopes. Such isotopes can appear 
during nuclear burning in stellar matter, particularly, in 
dense matter of white dwarf cores and neutron star en- 
velopes. The calculations were performed using the Sao 
Paulo (SP) potential in the context of the barrier pen- 
etration model; we employed the NL3-parametrization 
of nuclear density distributions of reactants within the 
Relativistic Hartree-Bogoliubov (RHB) approach. The 
model is based on the standard partial wave decomposi- 
tion (t — 0, 1, . . .) and considers motion of the nuclei in 
the effective potential 

h 2 f<£ + 1) 

U eS (r, E) = U c (r) + U SP (r : E) + \ (18) 

where Uc(r), Usp(r,E) and the last term arc the 
Coulomb, nuclear and centrifugal potentials, respectively. 
At low energies (E < Ec), the main contribution to S(E) 
comes from the £ = (s-wave) channel. The calcula- 
tional scheme is parameter-free and relatively simple for 
generating a set of data for many non-resonant reactions 
involving different isotopes. 

The reactions in question are summarized in Table [I] 
All isotopes studied were even-even nuclei. We consid- 
ered 10 reaction types, such as C+C and O+Ne, with the 
range of mass numbers for both species given in Table HI 
For each reaction, we computed S(E) on a dense grid 



TABLE I: Fusion reactions (Ai, Zi) + (A2, Z2) under consid- 
eration (after Ref. p|) 



Reaction 


Ai 


A 2 




Nr. of 


type 


even 


even 


MeV 


cases 


C+C 


10-24 


10-24 


17.9 


36 


c+o 


10-24 


12-28 


17.9 


72 


C+Ne 


10-24 


18-40 


19.9 


96 


C+Mg 


10-24 


20-46 


19.9 


112 


O+O 


12-28 


12-28 


19.9 


45 


O+Ne 


12-28 


18-40 


21.9 


108 


O+Mg 


12-28 


18-46 


21.9 


126 


Ne+Ne 


18-40 


18-40 


21.9 


78 


Ne+Mg 


18-40 


20-46 


24.9 


168 


Mg+Mg 


20-46 


20-46 


29.9 


105 



of E (with the energy step of 0.1 MeV) from 2 MeV to 
a maximum value i? max (also given in Table |l| covering 
wide energy ranges below and above the Coulomb bar- 
rier. The last column in Table U presents the number of 
considered reactions. 

The results of calculations using the SP model have 
been compared previously 0, 0, E3] with experimental 
data (if available) as well as with theoretical calculations 
performed using other models such as coupled-channels 
and fermionic molecular dynamics ones. As detailed in 
floj |. the calculated values of S(E) are uncertain due 
to nuclear physics effects - due to using the SP model 
with the NL3 nucleon density distribution. For the reac- 
tions involving stable nuclides, typical uncertainties are 
expected to be within a factor of 2, with maximum up 
to a factor of 4. For the reactions involving unstable nu- 
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TABLE II: Fit parameters of S(E) for reactions (Ai,Zi) + {A2,Z2) under consideration 



Reaction 


R 


AR la 


AR 2a 


ARib 


AR 2b 


8 


So 


Co 


£i 


Max. 


Rms 


type 


fm 


fm 


fm 


fm 


fm 




MeV b 






dev. 


dev. 


c+c 


7.4836 


0.1759 


0.1759 


0.0040 


0.0040 


0.0400 


1.3736 


3.5499 


0.2658 


0.35 


0.11 


c+o 


7.8671 


0.1740 


0.1280 


-0.0045 


-0.0310 


0.0412 


1.5438 


5.2576 


0.2306 


0.47 


0.12 


C+Ne 


7.9387 


0.1720 


0.1206 


-0.0171 


-0.0035 


0.0400 


1.9478 


3.7661 


0.2328 


0.62 


0.16 


C+Mg 


8.0513 


0.1705 


0.1014 


-0.0210 


-0.0186 


0.0386 


2.4327 


4.0059 


0.1844 


0.56 


0.17 


o+o 


8.0641 


0.1266 


0.1266 


-0.0377 


-0.0377 


0.0388 


2.1998 


6.0147 


0.1547 


0.50 


0.14 


O+Ne 


8.1191 


0.1257 


0.1183 


-0.0461 


-0.0068 


0.0371 


2.9486 


3.5127 


0.1702 


0.65 


0.19 


O+Mg 


8.2404 


0.1246 


0.0994 


-0.0500 


-0.0216 


0.0357 


3.7433 


2.8303 


0.1417 


0.65 


0.20 


Ne+Ne 


8.1419 


0.1175 


0.1175 


-0.0107 


-0.0107 


0.0348 


4.2215 


0.1225 


0.1717 


1.00 


0.27 


Ne+Mg 


8.2880 


0.1160 


0.0987 


-0.0157 


-0.0273 


0.0339 


5.2525 


0.2141 


0.1342 


0.86 


0.28 


Mg+Mg 


8.4509 


0.0976 


0.0976 


-0.0288 


-0.0288 


0.0332 


5.9785 


0.5263 


0.1393 


0.82 


0.28 



clei, typical uncertainties were roughly estimated to be 
as large as one order of magnitude, reaching two orders 
of magnitude at low energies for the reactions with very 
neutron-rich isotopes. These uncertainties reflect current 
state of art in our knowledge of S(E). 

B. Fits 

In Ref. @ we fitted the calculated S(E) by a 9- 
parameter phenomenological analytic expression. These 
fits are accurate (with maximum relative errors less than 
10%) but their use requires extensive tables (of 9 x 946 = 
8514 parameters). Here, we employ our new fit expres- 
sions (Sec. HI| and show that the same data can be ap- 
proximated using only 90 fit parameters. 

Let us consider reactions of each type (each line in Ta- 
ble H| separately and apply our analytic model (O and (|6]) 
to every reaction. In this way we determine 4 fit param- 
eters, 5*o, Ec, 5 and £, for every reaction. For instance, 
we have 4 x 105 = 420 parameters for Mg+Mg reactions. 
However, we notice that we can put So and 5 constant for 
all reactions of a given type (for instance, Sq — 5.9785 
MeV b and 6 = 0.0332 for all Mg+Mg reactions); this 
does not increase essentially the fit errors. Such constant 
So and 5 are given in Table [TTJ 

Still, we need to specify two parameters, Ec and £, 
for every reaction. Collecting the values of Ec and £ for 
all reactions of each type, we were able to fit them by 
analytic expressions 

E c =a/R { °\ (19) 

R { ° } = R + ARt \At - Aio| + Ai? 2 \A 2 - A 20 |; 
£ = £o + £i(Ai + A 2 ), (20) 

where A\q = 2Z± and A 2 q = 2Z 2 are mass numbers 
of most stable isotopes; ARi = AR la at A 1 > A w ; 
AR 1 = AR lb at At < A 10 ; AR 2 = AR 2a at A 2 > A 20 ; 
AR 2 = AR 2 b at A 2 < A 2 q. Thus, we have seven new fit 



parameters R, ARi a , AR 2a , ARu, Ai? 2 fc; £0 and £1 (also 
given in Table [n| for each reaction type, and, hence, 9 
parameters in total. Naturally, we have AR 2 = ARi for 
the reactions involving isotopes of the same clement (e.g., 
Mg+Mg). 

The eleventh column of Table [H] gives maximum rel- 
ative deviation of fitted S(E) from calculated ones for 
all reactions of a given type over all energy grid points 
(e.g., over 105 x 280 = 29400 points for the Mg+Mg re- 
actions). We see that the fitted values of S(E) do not 
deviate from the calculated ones by more than 100%. 
Root-mean square (rms) relative deviations, given in the 
last column, are a factor of 3-4 lower than the maxi- 
mum ones. This fit accuracy is acceptable because it 
is well within nuclear physics uncertainties of calculated 
S-factors fSec. HHAl) . 



C. Discussion 

Let us outline the main features of our fits for all 10 
reaction types. 

We start with the Mg+Mg reactions (Fig. [3]). They are 
characterized by the largest strength of Coulomb inter- 
action (largest product Z\Z 2 and largest a). On the left 
panel of Fig. [3] we compare calculated and fitted S(E) 
for six selected Mg+Mg reactions. Similar comparison 
for reactions of other types is given below. The solid 
curves are our fits; filled dots are calculated S(E). Here 
and below we plot the calculated data on a rarefied grid 
(with the energy step of 1 MeV) to simplify the figures. 
Open dots refer to E = Ec as determined from Eq. (| 19|) . 
For each reaction type we present six S(E) curves cho- 
sen in the following way. The lower curve corresponds 
to the reaction with the lightest isotopes ( 20 Mg+ 20 Mg 
in Fig. [3]). The upper curve is for the most massive iso- 
topes ( 46 Mg+ 46 Mg in Fig. GO). The second curve from 
the bottom is for the most stable isotopes ( 24 Mg+ 24 Mg); 
the third curve is for the lightest isotope 1 and the most 
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massive isotope 2 ( 20 Mg+ 46 Mg); and two next curves 
are for somewhat heavier isotope 1 and lighter isotope 2 
( 30 Mg+ 40 Mg; 40 Mg+ 40 Mg). The general trend is: the 
higher the reduced mass \x of the reacting nuclei, the 
larger S{E). We see that our 5(£')-model reproduces 
the data reasonably (uniformly) well for all Mg+Mg re- 
actions. All these fits arc done with the same 6 — 0.0332 
and Sq= 5. 9 78 5 MeV b. Fitting the same data with the 
model (fT6|) of rectangular potential would be much less 
accurate. This point is also illustrated on the left panel of 
Fig. [3] The dotted lines are our best fits of the same data 
with the rectangular potential. Such a potential leads to 
a faster growth of S(E) at E — > 0; it strongly (up to sev- 
eral orders of magnitude) overestimates low-energy S(E); 
fit errors become much higher. 

In the right panel of Fig. [3] we display the effective 
potential U(r) for the 46 Mg+ 46 Mg reaction (involving 
the most massive Mg isotopes from our collection). The 
solid line is our model U(r), reconstructed by fitting cal- 
culated values of S(E) with our model. It is given by 
Eq. and plotted also in Fig. [TJ with the fit parame- 
ters deduced from our fits (Table HI)) . The dotted line is 
our model U(r) for the best fit with 6 = 0. The long- 
dashed and short-dashed lines are the effective potentials 
U e g(r, E) given by Eq. (fl8|) and used in original SP cal- 
culations. They depend on E and are plotted for £7=10 
and 2 MeV, respectively. We see that the reconstructed 
potential with 8 — 0.0332 is remarkably close to the orig- 
inal ones. Therefore, by fitting the available S(E) data 
(calculated or experimental ones) with our 5(-E)-model, 
one can reconstruct the effective potential U(r). On the 
other hand, the model with S — gives the potential (the 
dotted line in the right panel of Fig. [3]) with unreasonably 
high Eq = 20.14 MeV that is sharply truncated at too 
large Rc = 10.29 fm. Naturally, this potential strongly 
overestimates S(E) at low E. 

The value of 5(0) for the 46 Mg+ 46 Mg reaction is ap- 
proximately 83 orders of magnitude higher than our fit 
value 5o=5.9785 MeV b. This huge difference is solely 
attributed to the definitions of 5(0) and Sq. While Sq is 
related to the penetrability of the real barrier (truncated 
at small r), 5(0) is defined through the penetrability of 
the barrier which remains pure Coulomb to r — > 0. These 
penetrabilities are drastically different. 

Figure 2] gives nine plots, similar to that on the left 
panel of Fig. [3j for the reactions of nine types. These 
reactions are ordered (from Ne+Mg, Ne+Ne and O+Mg 
in the top row to C+Ne, C+O and C+C in the bottom 
row) in such a way to have progressively lower strength of 
Coulomb nucleus- nucleus interaction (lower Z\Z-i). Gen- 
erally, the fits seem satisfactory and reasonably uniform. 

It is remarkable that the fit parameter Sq takes on the 
values in a narrow range from w 1.4 to ~ 6 MeV b for all 
reactions of our study, while the astrophysical 5-factor 
varies over many orders of magnitude. Our approach 
opens a possibility to extrapolate the values of Sq to a 
wider class of reactions without performing new calcu- 
lations of S(E). One can also extrapolate the values of 



S, Eq and £, and obtain thus 5(£7)-factors for new reac- 
tions. Note that we could have fitted all the data (Table 
HJ by using one and the same (5 « 0.04, and the fit accu- 
racy would be nearly the same as in our present fits. 

Our model S(E) is flexible to describe different S(E) 
curves. We believe that the description of S(E) in terms 
of four well defined parameters (Sq, 8, Ec, and £) is 
physical and sufficient for the majority of applications. 
However, our interpolation of Ec and £ by Eqs. (fT^|) 
and (|20p can be regarded only as a reasonably successful 
phenomenological fit. We expect that, while doing more 
accurate fitting of 5 (E) for these or other reactions in the 
future, one can find better (and physically meaningful) 
interpolation expressions for Ec and £ as functions of A\ 
and Ai, and consider S and 5o as functions of Ai and 
A2 as well. One can also improve our fit at E > Ec 
by going beyond the semi-classical approximation and 
by replacing the phenomenological S(E) dependence ([6]) 
with the results of more accurate consideration. We think 
that the 5(-E)-factor at E > Ec can be calculated in a 
more rigorous form and expressed through Ec, Sq and 8 
without introducing an additional parameter like £. 

Note that for each reaction type we consider some 
reactions involving proton-rich nuclei (e.g., 10 C+ 10 C) 
and many reactions involving neutron- rich nuclei (e.g., 
24 C+ 24 C). For each reaction type, we have observed a 
change in the behavior of fit parameters on Ai and A2 
while crossing the stability line (A ps 2Z). Because we 
include only a few proton-rich isotopes, we do not rec- 
ommend to extrapolate our fits to the region of proton- 
rich nuclei (that would require calculations of 5-factors 
for more proton-rich isotopes). We stress that our data 
sets include only even-even nuclei. In the future we can 
additionally calculate the 5-factors for reactions involv- 
ing other nuclei (even-odd or odd-odd) and approximate 
them with our model in a similar fashion. We do not rec- 
ommend directly extrapolating our present results (Table 
HT]) to these reactions. 



IV. 12 C+ 12 C REACTION 

Finally, let us discuss the quality of our S(E) model for 
the 12 C+ 12 C reaction. The reaction is very important for 
astrophysics of compact stars - for late stellar burning, 
ignition of type la supernovae and t rigg ering explosive 
events such as superbursts (e.g., [Til. Il2j] and references 
therein). Our results are summarized in Fig. [SJ On the 
left panel we show 5(£7)-data; on the right panel we give 
corresponding effective potentials U(r). 

The 12 C+ 12 C reaction cross section has been measured 
by many groups. The cross section contains resonances 
and the non-resonant part. We can compare our the- 
oretical model with the non-resonant contribution. We 
take experimental S(E) from Refs. @, 0, [HUH- The 
data cover the energy range from «2.1 MeV to 20 MeV. 
The data are nonuniform and not fully consistent, es- 
pecially at lower E, where experimental S(E) measure- 
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FIG. 3: (Color online) Left: 5*-factors for six Mg+Mg reactions. Filled dots are the Sao Paulo (SP) calculations; solid lines are 
our fits (Table [TTJ ; open dots show the fit values of Ec; dotted lines are best fits of S(E) assuming 5 = 0. Right: The effective 
potential U(r) for the 46 Mg+ 46 Mg reaction. The solid line, marked as Fit, is reconstructed from the calculated S(E) using our 
model (0 for U(r); the dotted line is the same but assuming rectangular potential (5 = 0). The long-dashed and short-dashed 
lines show the effective potential used in the original SP calculations of S(E) at E =10 and 2 MeV, respectively. 



merits are most difficult and experimental uncertainties 
are high. Experimental values of S(E) seem reasonably 
accurate at E > 4 MeV, moderately uncertain in the 
3 — 4 MeV range, and rather uncertain at lower E. The 
presence of low-energy resonances (e.g., 0, HH and ref- 
erences therein) complicates data analysis. To compare 
with the SP theoretical calculations, which neglect the 
resonances, and to smooth out the effect of experimen- 
tal uncertainties we bin the experimental data (with the 
bin width AE = 1 MeV) . The binned data are shown by 
crosses in the left panel of Fig. Theoretical SP S(E) 
values are presented by filled dots. 

The short-dashed line in the left panel of Fig. [5] is our 
4-parameter best fit of experimental S(E) factors. In this 
case we have E c = 6.25 MeV, 8 = 0.0763, S = 0.735 
MeV b, and x — 13.05. Reduced chi-square (per one 
degree of freedom) is 0.41, which means that the fit is 
acceptable. The solid line is our fit to the SP data. It is 
seen to be in reasonable agreement with the experimen- 
tal data. It agrees also with the fit to the experimental 
data at those energies at which the data are reliable. The 
dotted curve is another fit to the SP data, this time as- 
suming a rectangular barrier (1161) (5 = 0). It is in poor 
agreement with the experimental and SP data. 

Our 4-parameter fit to SP data gives (Table [TTJ) Ec — 
6.93 MeV, 5 = 0.04, S = 1.37 MeV b, x = 9.93. At 
subbarrier energies this fit is well described by Eqs. (|11[) 
and p3[) . with the expansion terms gi given by (|14p and 



(fP5). This yields 



S S p{E) = 1.7 x 10 lb exp(-0.302£: - 0.021^) MeV b, 

(21) 

where E is in MeV. Note that the contribution gu of 
the parabolic segment of the potential barrier U(r) to 
9i = 9ir + Su is generally substantial. For instance, 
go = 42.615 — 5.555 = 37.060. This expansion term 
is absorbed in 5(0) = 1.7 x 10 16 MeV b; the value 
goi = —5.555 reduces considerably 5(0). In the next 
term g 1 = -1.045 + 0.743 = -0.302 the parabolic U(r) 
part (gu = 0.743) is very essential. In the quadratic term 
g 2 = -0.023 + 0.002 = -0.021 it is less important, and 
the term itself is rather unimportant in (|21[) . 

A careful analysis of experimental data of different 
groups has recently been performed by Aguilera et al. 



21]. The authors rescaled some data trying to ob- 
tain a unified description of S(E) including resonant 
structures. Then they took several theoretical models 
of non-resonant S nv (E) and compared them with their 
full unified experimental S(E) data. The difference 
S(E) — S nT (E) was treated as the resonance contribu- 
tion. The best theoretical S ns (E) model was claimed to 
be that obtained with the KNS (Krappe-Nix-Sierk (23j ) 
barrier potential L^knsM- H gives most reasonable re- 
production of resonant structures after subtracting non- 
resonant term. 

It is customary (e.g., 0, |, 0, El [H| ) to approximate 
the non-resonant 5-factor for the 12 C+ 12 C reaction at 
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FIG. 4: (Color online) Each of nine panels shows S'-factors for six reactions of one type (as in the left panel of Fig. |3). The 
reaction types are Ne+Mg, Ne+Ne, O+Mg, O+Ne, C+Mg, O+O, C+Ne, C+O, and C+C. 



subbarrier energies by 

S m {E) = S a exp(-0.46£;), (22) 

where S a is a parameter [analogous to S, defined by 
(fl"2|) . with the specific value of g\ = —0.46 discussed 
below]. The results of Aguilera et al. [2l| can be ap- 
proximated in this way using their Fig. 11 [in our no- 
tations, that figure gives S(E) = S nr (E) exp(0.46£')]. 
Their best S nl (E) (given by the KNS model) in the en- 
ergy range from ss 3.5 to ~ 5.5 MeV (below the barrier 
but at those E where experimental data are reliable) can 
be approximated by (JH]) with S a « (1.4 - 1.7) x 10 16 
MeV b (earlier studies gave S a «3x 10 16 MeV b; e.g., 



24]). However, this approximation becomes inaccu- 
rate at lower E. We remark, that at any E < 5.5 MeV, 
the best (KNS) model of Ref. [23| is accurately approx- 
imated by Skns(E) w 0.83 x 10 16 exp(-0.32£:) MeV b, 
in reasonable agreement with our approximation (|21[) of 
SP data. Nevertheless, another acceptable S nr (E) model 
of Ref. 1211 . based on the proximity-adiabatic (PA) ap- 
proach [25], is approximated at the same energies E < 5.5 
MeV as S PA (E) « 1.6 x 10 16 exp(-0.46E) MeV b, in 
agreement with (|22|) . 

These different approximations reflect uncertainties in 
our knowledge of S m (E). One should be careful in us- 
ing for the 12 C+ 12 C reaction. The factor g x in the 
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FIG. 5: (Color online) Left: S'-factors for the 12 C+ 12 C reaction. Crosses show experimental data, short-dashed line is the 
best-fit. Filled dots are SP calculations, solid line is our fit with 5 — 0.04 (Table [TTJ ; dotted line is our best fit with 8 — 0. 
Right: The effective potential U(r) for the same reaction. Short-dashed line is reconstructed from the fit to experimental data; 
long-dashed line is the original SP potential; solid and dotted lines are the potentials U(r) reconstructed from the fits to SP 
data with 5 — 0.04 and 0, respectively. See text for details. 



exponent argument can be different from 171 = —0.46 
(can be closer to —0.3). The value g\ = —0.46 was first 
introduced by Patterson et al. 0- The authors claimed 
that it came from the model of rectangular potential bar- 
rier. According to Eq. (1171) . in this model we would have 
3/2 

<7irct = —0.05 Rfj , where Rc = Rci is m fm- Taking 
the standard value Rc ~ 2 x 1.3 A 1 / 3 = 6 fm we would 
get <7i rc t = —0.73, noticeably different from —0.46. Tak- 
ing gxrct — —0.46 one would have Rc — 4.4 fm [6j, an 
unrealistically small radius of the potential well U(r). In 
fact, Patterson et al. [6j obtained g\ = —0.46 by fitting 
a restricted set of experimental S(E) values available by 
1969. They discussed possible variations of g\ but the 
discussion has been forgotten, whereas their best value 
gi = —0.46 is widely cited in the literature. 

The right panel of Fig. [5] presents the effective po- 
tentials U(r) corresponding to S(E) models on the left 
panel. The short-dashed curve gives U(r) that is cal- 
culated from Eq. ([7]) with the parameters Ec = 6.25 
MeV and 5 = 0.0763 inferred from the fit to experimental 
S(E). It is our reconstruction of the real potential U(r) 
(whose details are still unknown) using our analytic S(E) 
model. The long-dashed line is the theoretical SP U e g(r) 
(that is almost independent of E for the 12 C-r- 12 C reac- 
tion); it was used to calculate the SP S(E). The solid 
curve is U(r) reconstructed from the full (Ec = 6.93 
MeV, S = 0.04) fit to our calculated SP S(E). The dot- 
ted curve is a similar reconstruction but based on the 
rectangular U(r) model [Eq. (fTo) . (5 = 0]. The dotted 
curve looks unrealistic - it gives too large Vc — 9.18 
MeV, indicating once more that the model of rectangu- 



lar potential is too crude. Three other potentials are re- 
markably close even in this, most difficult 12 C+ 12 C case, 
complicated by pronounced resonance structures of ex- 
perimental S(E). This fact confirms that the SP model 
is generally a valid tool for studying non-resonant fusion 
reactions. Moreover, we see that our analytic S(E) model 
can help to reconstruct U (r) from experimental data. Let 
us add that the solid, short-dashed and long-dashed U(r) 
curves are also close to effective potentials, particularly, 
to C/KNsMand C/paO"), used by Aguliera et al. [2l[ to 
approximate non-resonant contribution to S(E) in ex- 
perimental data. As seen from Fig. 8 and Table 3 in [2l[ , 
the basic parameters (Ec and Rc) of t/KNs(V)and J/paW 
are fairly close, but Ukns ( r ) has a slightly more extended 
low-r wing which, however, changes S nT (E) behavior at 
low E (from <7ikns = —0.32 to (?ipa = —0.46). 

Clearly, different segments of S(E) are determined by 
different parts of U(r). The range of E below Ec down 
to 3 — 4 MeV is controlled by U(r) at r > 6 fm [not far 
from the U(r)-peak}. In this case, different theoretical 
J7(r)-models give reasonably similar non-resonant S(E)- 
factors, which generally agree with (accurate) experimen- 
tal data. At lower energies, E ~ 1 — 3 MeV, which are 
important for 12 C burning in stellar matter in the ther- 
monuclear regime, experimental values of S(E) are either 
uncertain or not available. These values are sensitive to 
the sharp low-r wing of the U(r) potential (r ~ 5 — 6 
fm) which is not very well constrained by theory and ex- 
periment. However, because the 1 — 3 MeV energy range 
is close to the range of higher E, where the S(i?)-factor 
is well studied, one expects that an extrapolation from 
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higher E to the 1 — 3 MeV range is more or less reli- 
able. Finally, lowest energies E < 1 MeV are important 
for pycnonuclear burning of 12 C in dense stellar matter 
(e.g., [4, 9]). In this case, S(E) is controlled by the very 
steep low-r slope of U(r) and seems rather uncertain. It 
can be affected by slight variations of the U(r) slope (as 
discussed above taking KNS and AP models in [2l| as 
an example). Extrapolations to these energies can be 
inaccurate. The problem is further complicated by the 
oblateness of 12 C nuclei in the ground state. Our analysis 
is based on the approximation of spherically symmetric 
nuclei and spherically symmetric potential U(r). In the 
presence of oblateness, 5(.E)-factors depend on orienta- 
tions of colliding nuclei. This effect is beyond the scope of 
the present paper. Note, however, that according to cal- 
culations [26| the oblateness increases low-energy S(E) 
for the 12 C+ 12 C reaction by a factor of ~ 1.7 [which is 
within theoretical uncertainties of non-resonant S(E)]. 

We have also compared calculated and fitted S(E)- 
factors with experimental data for the 12 C+ 16 and 
16 0+ 16 reactions. We have carried out SP calculations 
based on the RHB approach outlined in Sec. MI Al as well 
as SP calculations which employ two-parameter Fermi 
(2pF) parametrization of nuclear density distributions of 
reactants. 2pF calculations agree with experimental data 
better than RHB ones (although the accuracy of both 
approaches is sufficient for many applications). If we re- 
stricted ourselves to reactions involving stable nuclei, the 
2pF parametrization would be more accurate. However, 
our main goal was to obtain a uniform set of theoretical 
S(E) data for a large collection of nuclei involving un- 
stable ones (to simulate nuclear burning in neutron stars 
and white dwarfs). In this case, SP calculations based on 
the RHB approach are favorable. 

V. CONCLUSIONS 

We have suggested (Sec. [TTJ) a simple model with phys- 
ically meaningful parameters to describe the astrophysi- 
cal 5-factor as a function of center-of-mass energy E of 
reacting heavy nuclei for non-resonant fusion reactions. 
Our main conclusions are as follows: 

• For any reaction, the model gives S(E) in an an- 
alytic form in terms of four parameters. They are 
Ec, the height of the Coulomb barrier; Sq that 
characterizes the strength of the nucleus-nucleus 
interaction neglecting Coulomb interaction; d that 
describes the peak broadening of the effective bar- 
rier potential U(r); and £ to describe the transition 
from subbarrier energies to E 3> Ec- The model 
is expected to be sufficiently accurate for energies 
below and above Ec (up to a few Ec)- 

• As an example, we have applied our model to de- 
scribe the 5-factors for 946 fusion reactions involv- 
ing various isotopes of C, O, Ne, and Mg, from 
the stability valley to very neutron-rich nuclei. In 



Ref. [|| these ^-factors were calculated using the SP 
method and the barrier penetration model. They 
were fitted by a phenomenological formula con- 
taining 9 fit parameters for every reaction (8514 
parameters in total). With the present analytic 
model, we can fit the same data set using 90 fit 
parameters fSec. IIIIBl Table HI]) . The fit accuracy 
is worse than in Ref. [5j but is well within esti- 
mated nuclear-physics uncertainties of calculated 
S{E) (Sec. EE). 

• We have also compared our model (Sec. IIVI) for 
the 12 C+ 12 C reaction (that is most important for 
neutron stars and white dwarfs) with experimental 
S(E) data and discussed the problem of extrapola- 
tion of experimental S(E) to low energies of astro- 
physical importance. 

• Our analytic S'(£ , )-model is easy for implementing 
into computer codes, which calculate nuclear re- 
action rates and simulate various nuclear burning 
phenomena in astrophysical environment; it is not 
as costly for CPU time as reading large tables can 
be. The S(E) dependence for any reaction is de- 
termined by the values of Ec , So , 5 and £ through 
Eqs. ([3]), ^ and ([§]). As a rule, one needs only 
subbarrier S(E) to calculate nuclear reaction rates 
in stellar matter. In this case, it is sufficient to use 
a simplified expression (|11[) . where the coefficients 
gi, <?2j and 53 are given by Eqs. (fT4)) and (fT5j). 

• The analytic model is practical for describing large 
uniform sets of S(E) data (for instance, many reac- 
tions involving isotopes of the same elements) . The 
parameters Ec, Sq, S and £ vary slowly from one 
reaction to another, and are easily interpolated (by 
analytic expressions) over large sets; So and 5 can 
be set constant for many reactions. Analytic inter- 
polations can be used to extrapolate S(E) to other 
reactions of the same type. 

• The functional form of our analytical S(E) is flex- 
ible to describe qualitatively different behaviors of 
S(E). Particularly, by varying d we can obtain ei- 
ther growth or decrease of S(E) as E — > 0. The 
decrease is realized at not too small 5 and may ex- 
plain the low-energy hindrance of S(E) whose sig- 
nature was observed in some reactions (e.g., [Tll.[27| 
and references therein). The low-energy behavior 
of S(E) is indeed very sensitive to the parameter 5. 
The model of sharply truncated Coulomb poten- 
tial (fP6|) . that is widely used in the astrophysical 
literature (e.g., Refs. 0, [|j]), can be inaccurate in 
extrapolating calculated or measured S(E) to low 
E (Sees. HHC] and EJ. 

• Fitting a given S (E) (computed or measured in lab- 
oratory) with our analytic model can be used to re- 
construct the effective potential U(r) (Figs. |3] and 
[5]). Of course, the real potential U(r) can be too 
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complicated to be exactly described by our model 
potential 0. However, we expect that this poten- 
tial allows one to reproduce correct S(E) behavior 
in many cases. 

There is no doubt that the Z7(r)-peak is not sharp, 
but broadened. Roughly speaking, this broadening is 
twofold. First, the potential peak becomes smoother. 
Second, the low-r wing of U(r) becomes less steep. 
Clearly, the second effect has stronger impact on the low- 
energy S(E) than the first one. In our model, both effects 
are described by one and the same parameter 5. We can 
complicate the model by introducing new parameters but 
think that the present version is good as the first step. 
Let us add that our model is useful for reactions between 
heavy nuclei. Astrophysical S'-factors for reactions in- 
volving light nuclei contain strong resonances (e.g., [241 ]) 
which are not described by our model. 

Let us add that nuclear reaction rates in dense stellar 
matter (especially, in the cores of white dwarfs and en- 
velopes of neutron stars 28]) can be greatly affected by 
plasma screening of the Coulomb interaction and by the 
transition to pycnonuclear burning regime (where zero- 
point vibrations of nuclei in a strongly coupled plasma 
of ions become important). These plasma physics effects 
were described by Salpeter and Van Horn [29j (also see 
0) B Hfll, HU and references therein). They modify the 
interaction potential U(r) but mainly at sufficiently large 
r, typically higher than nucleus sizes, while we focus on 
the nuclear physics effects which influence U(r) at lower 
r. It is widely thought that the plasma physics and nu- 
clear physics effects arc distinctly different and can be 
considered separately. However, we notice that in very 



dense and not too hot stellar matter both effects can be- 
come interrelated (and should be considered together). 

We expect that the broadening of the U(r) peak is es- 
pecially important for pycnonuclear reactions in the in- 
ner crust of accreting neutron stars in X-ray transients 
[32M34j ]. They are compact binary systems containing a 
neutron star and a low-mass companion. Pycnonuclear 
reactions are thought to be responsible for deep crustal 
heating of accreted matter. The heating can power (35[ 
thermal surface emission of these neutron stars that is 
observed in quiescent states of transients (see, e.g., Refs. 
[ID, [36[). Pycnonuclear reactions occur at high densities 
and involve very neutron-rich nuclei (e.g;, 34 Ne+ 34 Ne at 
p w 2 x 10 12 g cm" 3 , according to Ref. [33[) immersed in 
a sea of free neutrons (e.g., [33]). The [/(r)-peak should 
be broadened not only by a diffusive structure of neutron- 
rich nuclei (that is taken into account in the SP calcu- 
lations) but also by the presence of free neutrons. The 
latter effect is unexplored, but it can affect S(E), nuclear 
reaction rates, the deep crustal heating and its observa- 
tional manifestations (e.g., Ref. [38jV 
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